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HIGHLIGHTS 


•  Thermoelectric  SiC  nose-tip  structure  for  aerodynamic  heat  harvesting  of  high-speed  vehicles  is  studied. 

•  Thermoelectric  performance  is  predicted  based  on  numerical  methods  and  experimental  thermoelectric  parameters. 

•  The  effects  of  Mach  number  on  thermoelectric  performance  are  studied  in  the  present  paper. 

•  Results  with  respect  to  the  Thomson  effect  are  also  explored. 

•  Output  power  and  energy  efficiency  of  the  thermoelectric  nose-tip  are  increased  with  the  increase  of  Mach  number. 
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This  paper  focus  on  the  effects  of  Mach  number  on  thermoelectric  energy  conversion  for  the  limitation  of 
aero-heating  and  the  feasibility  of  energy  harvesting  on  supersonic  vehicles.  A  model  of  nose-tip 
structure  constructed  with  SiC  ceramics  is  developed  to  numerically  study  the  thermoelectric  perfor¬ 
mance  in  a  supersonic  flow  field  by  employing  the  computational  fluid  dynamics  and  the  thermal 
conduction  theory.  Results  are  given  in  the  cases  of  different  Mach  numbers.  Moreover,  the  thermo¬ 
electric  performance  in  each  case  is  predicted  with  and  without  Thomson  heat,  respectively.  Due  to  the 
increase  of  Mach  number,  both  the  temperature  difference  and  the  conductive  heat  flux  between  the  hot 
side  and  the  cold  side  of  nose  tip  are  increased.  This  results  in  the  growth  of  the  thermoelectric  power 
generated  and  the  energy  conversion  efficiency.  With  respect  to  the  Thomson  effect,  over  50%  of  total 
power  generated  converts  to  Thomson  heat,  which  greatly  reduces  the  thermoelectric  power  and  effi¬ 
ciency.  However,  whether  the  Thomson  effect  is  considered  or  not,  with  the  Mach  number  increasing 
from  2.5  to  4.5,  the  thermoelectric  performance  can  be  effectively  improved. 

©  2013  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

When  cruising  in  atmosphere,  supersonic  vehicles  are  con¬ 
fronted  with  high  thermal  loads  on  leading-edges.  There  have  been 
many  efficient  heat  dissipation  mechanisms  for  preventing 
extremely  dangerous  damage  to  the  geometry  configuration  of 
vehicles  from  serious  aero-heating,  such  as  blunted  waverider  [1], 
film  cooling  [2],  heat-pipe  cooling  [3],  transpiration  cooling  [4]  and 
opposing  jet  cooling  [5].  However,  the  dispersed  aerodynamic  heat 
is  hardly  reusable  and  recycled. 
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Recently,  a  feasible  application  was  proposed  for  the  supersonic 
vehicles  aiming  at  converting  aero-heating  quantity  to  electricity 
by  thermal  photovoltaic  system  (TPV)  [6].  TPV  is  a  kind  of  ther¬ 
moelectric  energy  conversion  technique,  which  enables  heat- 
electricity  direct  conversion.  Instead  of  dispersing  most  of  aero¬ 
dynamic  heat,  the  thermoelectric  technology  enables  the  energy 
harvesting  from  it.  Thermoelectric  technology  is  very  promising  in 
low  level  thermal  energy  recovery.  Jang  and  Tsai  [7]  studied  the 
optimization  of  thermoelectric  modules  attached  to  the  chimney 
plate  for  waste  gas  heat  recovery.  Kim  et  al.  [8]  designed  a  new 
thermoelectric  generation  system  with  enlarged  hot  surface  area 
for  heat  pipes.  Nuwayhid  et  al.  [9]  used  commercial  available 
thermoelectric  modules  to  evaluate  the  feasibility  of  stove-run 
thermoelectric  generators  for  domestic  electricity  supply.  Min 
and  Rowe  [10]  developed  a  thermoelectric  system  with  combusted 
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Nomenclature 

V 

thermoelectric  voltage  (mV) 

x,  y,  z 

cartesian  coordinates 

Alphabetization 

Greek 

symbols 

A 

area  (cm2) 

a 

Seebeck  coefficient  (V/K) 

Fj,  Gj 

inviscid  terms  for  flow 

V 

energy  conversion  efficiency  (%) 

EV,  FVl  Gy 

viscous  terms  for  flow 

K 

thermal  conductivity  (W/(m  I<)) 

I 

current  (A) 

TC 

Peltier  coefficient  (m) 

M 

Mach  number 

P 

electric  resistivity  (Q  mm) 

P 

output  power  (W) 

0 

heat  flux  rate  (W) 

Pt  h 

power  excluding  Thomson  heat  (W) 

P 

pressure  (kPa) 

Subscripts 

Q 

heat  transfer  rate  (W) 

C 

cold  side 

Qjh 

Thomson  heat  rate  (W) 

H 

hot  side 

R 

electric  resistance  (Q) 

L 

Load 

S 

source  terms  for  flow 

00 

free  stream  static 

St 

Stanton  number 

T 

absolute  temperature  (I<) 

Abbreviations 

AT 

temperature  difference  (I<) 

N-S 

Navier— Stokes 

t 

time  (s) 

SST 

shear  stress  transport 

U 

conservation  variable 

heat  as  its  thermal  input.  It  is  widely  recognized  that  the  thermo¬ 
electric  equipment  and  devices  are  reliable,  safe  and  capable  to 
generate  a  minimum  acceptable  power,  and  have  been  found  very 
extensive  applications  in  specific  areas,  such  as  military  and  aero¬ 
space  [11  ].  As  sunshine  and  other  light  sources  most  needed  for  the 
photovoltaic  units  may  not  be  available  in  the  cabin,  the  use  of  TPV 
may  not  be  as  efficient  as  expected.  Since  the  main  feature  of  the 
temperature  distribution  on  supersonic  vehicles  leading-edges  is 
the  large  temperature  difference  rather  than  the  multiple  thermal 
sources  provided,  it  is  more  feasible  to  use  structural  materials  with 
thermoelectric  function  [12]. 

SiC  matrix  ceramics  [13]  and  SiC-based  composites  [14]  are 
high-temperature  structural  materials  widely  applied  in  space  ve¬ 
hicles.  Besides,  SiC  ceramics  perform  good  thermoelectric  proper¬ 
ties  at  high  temperatures  [15],  and  can  be  applied  to  chip  cooling 
[16].  Even  though  the  thermoelectric  figure  of  merits  of  SiC  ce¬ 
ramics  is  relatively  low,  it  can  be  greatly  improved  by  several  ap¬ 
proaches,  for  example,  doping  [17].  So  far,  the  prospect  of  SiC 
materials  in  thermoelectric  applications  is  notable  for  their  thermal 
stability  and  high  specific  strength  in  large  temperature  range.  In 
this  work,  a  nose-tip  model  constructed  with  SiC  ceramics  is 
developed  to  predict  the  thermoelectric  energy  conversion  per¬ 
formance  in  a  supersonic  flow  field  by  computational  fluid  dy¬ 
namics  and  thermal  conduction  theory.  The  effect  of  Mach  number 
on  the  thermoelectric  performance  is  discussed  in  detailed. 

2.  Materials  and  methods 

2.1.  Nose-tip  model  and  numerical  methods 

As  shown  in  Fig.  1,  the  nose-tip  model  is  40.6  mm  in  length,  and 
has  a  hemisphere-shaped  leading  edge,  which  is  25.5  mm  in 
diameter.  The  flow  field  parameters  and  heat  fluxes  were  obtained 
by  solving  three-dimensional  N-S  governing  equations  with  the 
SST  turbulence  model  [18,19].  The  steady-state  full  N-S  equations 
in  conservation  form  are  given  by 

dU  dE[  6Fj  6Gi  d Ev  d Fv  dGv  .  . 

-w+^+^+ni-i»+w+^+s  m 

Because  no  chemical  reactions  are  involved,  S  is  set  to  zero.  The 
convective  terms  are  approximated  using  the  Van-leer  flux  splitting 


method  [20],  and  the  viscous  terms  using  the  central  difference 
method.  A  semi-explicit  Runge-Kutta  scheme  is  used  for  the  time 
integration  [21  ].  The  temperature  distributions  of  the  flow  field  and 
solid  nose-tip  model  are  obtained  by  solving  heat  conduction 
equation  [22].  The  grids  of  the  simulation  model  in  the  symmetrical 
plane,  outlet  and  the  surface  of  the  nose-tip  are  shown  in  Fig.  2.  The 
boundary  conditions  are  listed  in  Table  1.  It  is  assumed  that  the 
fluid-solid  interface  is  no-slip  and  the  nose-tip  solid  region  is 
unexpanded  and  undistorted.  The  radius  heat  transfer  in  solid  re¬ 
gion  is  ignored,  because  the  maximum  temperature  is  nearly 
1000  K,  below  which  the  radius  heat  is  relatively  small. 

2.2.  Evaluation  of  thermoelectric  performance 

In  an  effective  circuit  shown  in  Fig.  3,  the  nose-tip  structure 
serves  as  the  thermoelectric  module.  Only  the  thermoelectric 
properties  of  the  nose-tip  materials  (SiC)  are  considered  to  study 
specially  the  recovery  capability  of  model.  When  temperature  dif¬ 
ference  applying  on  the  nose-tip,  the  generated  open  circuit  voltage 
is  given  by 

V  =  nH-nc  (2) 

Assuming  constant  thermal  conductivity  and  internal  electric 
resistance  of  nose-tip,  the  rates  of  supply  heat  flux  at  the  hot  side 
and  removal  heat  flux  at  the  cold  side  are  expressed  as 

=  tch/  -  Q  -  0.5 12R  (3) 
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Fig.  2.  Grids  on  the  simulation  model. 


Fig.  3.  Schematic  diagram  of  nose-tip  thermoelectric  module. 


<f>c  =  7tc/  -  Q  +  0.5 12R  (4) 

of  which  the  items  on  the  right  side  in  turn  are  Peltier  heat, 
conductive  heat,  and  Joule  heat.  Instead  of  the  actual  value,  the 
average  internal  electric  resistance  is  expressed  as 

Th 

%cJh)  =  J  R(T)dT/(TH  -  7c)  (5) 

Tc 

The  generated  power  can  be  calculated  by  either  multiplying 
output  current  with  loaded  electric  resistance  or  the  variation  of 
heat  flux  as  follows: 

P  =  /2  rl  =  V2Rl/(R  +  Rl )2  =  <PH-<PC  (6) 

The  thermoelectric  energy  conversion  efficiency  of  nose-tip 
model  is  given  by 

V  =  P/$  H  (7) 

Thomson  heat  relates  to  the  rate  of  generation  of  reversible  heat 
with  results  from  the  passage  of  a  current  along  a  portion  of  a  single 
conductor  along  which  there  is  a  temperature  difference  AT. 
Although  the  Thomson  heat  is  not  of  primary  importance  in  ther¬ 
moelectric  devices,  it  must  be  considered  in  detailed  calculations 
[23,24].  Thomson  heat  and  the  output  power  are  illustrated  by 


Qjh  =  (duH/d T  -  dirc/dr)/AT 

(B) 

^Th  =  (UH  “  KC)I  ~  Qrh  “  !2R 

(9) 

Accordingly,  the  energy  efficiency  considering  Thomson  effect  is 
expressed  as 


Table  1 

Computational  boundary  conditions  of  flow  field  and  nose  region. 


Item  Value 


Free  stream  Mach  number,  M  2.5— 4.5 

Free  stream  static  pressure,  p*,  5.475  kPa 

Free  stream  static  temperature,  T x  217  K 

Cold  side  temperature  of  prototype,  Tc  300  K 


^Th  =  Pjh/^H  (10) 

2.3.  Thermoelectric  properties  of  the  SiC  materials 

Series  of  SiC  ceramics  and  composites  have  been  investigated  on 
their  thermoelectric  properties.  However,  there  are  few  compre¬ 
hensive  reports  on  the  thermoelectric  parameters  of  polymer- 
derived  SiC  materials  applied  for  aerospace  structures,  such  as 
thermal  conductivity,  electric  resistivity  and  Seebeck  coefficient.  In 
order  to  obtain  the  needed  parameters,  dense  SiC  ceramics  were 
prepared  with  the  preceramic  precursor  polycarbosilane  produced 
by  our  laboratory  [25].  The  Peltier  coefficient  and  electric  conduc¬ 
tivity  were  measured  by  four-probe  method  (ZEM-3,  ULVAC-RIKO, 
Japan).  The  thermal  conductivity  was  studied  using  laser-flashed 
method  (LFA457,  Netzsch,  Germany).  Archimedes  method  was 
used  to  measure  the  density  of  the  consolidated  SiC  sample  with 
deionized  water  as  the  immersion  medium. 

The  measured  density  of  the  consolidated  SiC  sample  is  1.72  g / 
cm3.  Since  the  thermal  conductivity  varies  very  slightly  in  the 
temperature  range  of  300  K-923  K,  it  is  assumed  to  be  constant, 
and  its  average  value  is  5.73  W/(m  K).  Peltier  coefficient  (tt)  and 
electric  resistivity  (p)  versus  temperatures  of  the  SiC  ceramics  are 


Fig.  4.  Peltier  coefficient  and  electric  resistivity  of  the  dense  SiC  ceramic  materials. 
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Fig.  5.  Experimental  and  numerical  Stanton  number  along  nose-tip  surface 
comparison. 


shown  in  Fig.  4.  With  the  elevated  temperature,  the  absolute  value 
of  tz  increases,  and  inversely,  p  decreases.  Results  suggest  good 
thermoelectric  properties  of  SiC  at  even  higher  temperatures. 

3.  Validation  of  computational  methods 

Prior  to  computation,  a  verification  of  the  algorithm  and  grid 
independence  of  the  numerical  solution  is  performed  to  ensure  the 
accuracy  and  validity  of  the  results.  Flow  conditions  is  chosen  ac¬ 
cording  to  the  related  experiment  in  Ref.  [26].  As  shown  in  Fig.  5, 
the  predicted  Stanton  numbers  (St)  along  the  surface  of  the 
hemisphere  nose  are  compared  with  the  experimental  data.  The 
numerical  and  experimental  results  are  in  good  agreement  with 
each  other.  There  are  nine  experimental  data  points,  except  the  first 
one,  the  simulation  accuracies  of  the  remaining  points  are  in  15%,  of 
which  the  best  one  is  close  to  1%.  The  discrepancy  is  attributed  to 
the  computation  error,  assumption  of  the  simulation  model  and 
experimental  measurement  error.  Especially  with  the  stagnation 
point,  the  discrepancy  is  much  bigger  for  two  main  reasons.  For  one 
thing,  due  to  a  significant  temperature  gradient  around  the  stag¬ 
nation,  a  little  change  of  position  will  cause  a  large  variety  of  the 
heat  flux;  for  the  other,  it  is  hard  to  measure  the  heat  flux  accurately 
in  the  stagnation  region. 

4.  Results  and  discussion 

By  means  of  above  numerical  procedures,  the  flow  parameters 
and  the  thermal  conditions  can  be  obtained.  The  Mach  number 
distribution  in  flow  field  (M  3.5)  is  shown  in  Fig.  6.  And  the  tem¬ 
perature  distributions  of  the  flow  field  and  nose  tip  are  given  in 
Figs.  7  and  8,  respectively.  The  highest  temperature  of  the  super¬ 
sonic  flow  is  over  746  K,  and  the  temperature  at  the  stagnation  of 
the  nose-tip  region  is  up  to  655  K.  With  the  temperature  at  the  cold 
side  of  the  nose  model  set  to  300  K,  as  a  result,  a  temperature 
difference  between  the  two  sides  of  355  K  is  obtained. 

Both  the  temperature  difference  AT  between  the  hot  and  the 
cold  surfaces  and  the  conductive  heat  flux  Q.  increase  as  the  Mach 
number  increases.  As  shown  in  Fig.  9,  at  M  4.5,  the  temperature 
difference  reaches  as  high  as  670  K,  nearly  seven  times  of  that  at  M 


Mach  Number 
M  3.  5 

I™ 

2 . 626e+000 

I 

i  1 .751e+000 

8 . 754e-001 
^ 1 . 000e-015 


Fig.  6.  Mach  number  distribution  in  the  flow  field  (symmetrical  plane,  M  3.5). 


2.5.  And  the  heat  flux  density  is  nearly  four  times  of  that  at  M  2.5. 
On  the  other  hand,  the  internal  electric  resistance  decreases,  while 
the  thermoelectric  voltage  of  the  nose-tip  model  increases  with  the 
increase  of  Mach  number,  as  shown  in  Fig.  10.  The  thermoelectric 
voltage  is  about  0.175  mV  at  M  4.5,  and  the  electric  resistance  of  the 
nose  model  reduces  to  nearly  0.12  Q.  Results  indicate  better  ther¬ 
moelectric  performance  at  high  Mach  numbers. 

The  performance  of  the  thermoelectric  nose  tip  can  be  charac¬ 
terized.  Under  the  conditions  of  different  Mach  numbers,  the 
generated  output  power  and  the  energy  efficiency  as  a  function  of 
current  neglecting  Thomson  effect  are  illustrated  in  Fig.  11.  From  M 
2.5  to  M  4.5,  the  maximum  current  increases  from  0.12  A  to  1.44  A. 
Larger  current  is  permitted  at  higher  Mach  numbers,  in  which  case 
the  circuit  has  larger  thermoelectric  voltage  and  smaller  internal 
electric  resistance.  When  the  loaded  electric  resistance  of  the  cir¬ 
cuit  equals  to  the  internal  electric  resistance  of  the  nose-tip,  the 
maximum  power  and  maximum  efficiency  is  obtained.  As  shown  in 


Temperature  (K) 
M  3.  5  flow  field 
7 . 461e+002 


6. 137e+002 
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Fig.  7.  Temperature  distribution  in  the  flow  field  (symmetrical  plane,  M  3.5). 
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M  3.  5  solid  field 
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Fig.  8.  Temperature  distribution  in  the  nose- tip  region  (symmetrical  plane,  M  3.5). 


Fig.  10.  Thermoelectric  voltage  and  internal  electric  resistance  of  nose-tip  at  different 
Mach  numbers. 


Fig.  12,  both  the  maximum  of  output  power  and  energy  efficiency 
are  increased  with  the  elevated  Mach  numbers  attributed  to 
enlarged  temperature  difference,  heat  flux,  thermoelectric  voltage, 
and  the  decrease  of  internal  electric  resistance  (also  see  Figs.  9  and 
10).  The  maximum  of  power  output  at  M  4.5  is  over  0.062  W,  30 
times  larger  than  that  of  M  2.5.  That  is  almost  0.25  W/cm2  on 
average  (divided  by  the  cold  side  area,  about  0.204  cm2).  Corre¬ 
spondingly,  the  maximum  energy  efficiency  increased  with  the 
increase  of  Mach  number.  At  M  4.5  the  maximum  efficiency  is 
0.22%,  about  10  times  of  that  at  M  2.5. 

Thomson  heat  would  grow  much  faster  if  there  were  large 
temperature  difference  or  exceeding  circuit  current.  Under  the 
condition  of  each  Mach  number,  the  output  power  reduces  to  less 
than  a  half  of  the  total  generated  power.  As  shown  in  Fig.  13,  at  M 
4.5  the  output  power  is  0.021  W,  which  possesses  less  than  30%  of 
the  generated  thermoelectric  power.  The  energy  conversion  effi¬ 
ciency  also  decreases  dramatically.  The  maximum  value  reduces  to 
0.0085%  at  M  4.5,  less  than  a  half  of  the  one  without  respect  to  the 
Thomson  effect  (0.022%).  It  is  also  observed  that  the  maximum 
current  is  reduced  with  respect  to  the  Thomson  effect.  This  is 
because  partial  generated  power  transfers  to  Thomson  heat,  and 
relatively  smaller  electric  power  is  remained. 

Although  higher  Mach  number  results  in  much  greater 
Thomson  heat,  the  thermoelectric  performance  will  be  improved 
due  to  the  continuously  increasing  thermoelectric  voltage  and 
decreasing  electric  resistance  of  the  nose-tip.  In  order  to  maintain 
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thermoelectric  voltage  large  enough,  the  absolute  value  of  See- 
beck  coefficient  (a)  should  be  enlarged.  On  the  other  hand,  ac¬ 
cording  to  (8),  of  which  the  Peltier  coefficient  is  related  to 
Seebeck  coefficient  a  by  Kelvin  relations: 

a  =  n /T  (11) 

the  Seebeck  coefficient  should  remain  constant  in  a  wide  range  of 
temperature  in  order  to  eliminate  the  effect  of  Thomson  heat.  It  is 
well  recognized  that  nanostructured  materials,  such  as  super¬ 
lattices,  quantum  dots,  nanowires,  and  nanocomposites  play  a  role 
in  optimizing  the  Seebeck  coefficient  [27].  For  this  reason,  the 
application  of  the  thermoelectric  nose-tip  for  aerospace  vehicles  at 
high  Mach  numbers  is  encouraging. 

It  is  notable  that  the  thermal  conductivity  and  electric  resistance 
are  assumed  constant  at  each  case,  and  that  the  thermoelectric 
generated  model  is  reduced  to  one  dimensional  problem.  The 
computational  results  may  not  be  in  good  agreement  with  the  real 
ones.  With  further  studies  the  three-dimensional  numerical  model 


Fig.  9.  Temperature  difference  between  the  hot  and  cold  side  and  conductive  heat  flux  Fig.  11.  Output  power  and  energy  efficiency  of  nose-tip  versus  current  at  different 
density  at  different  Mach  numbers.  Mach  numbers. 
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introduced,  the  results  will  be  more  accurate.  Moreover,  the 
collection  of  experimental  data  is  a  need  for  the  practical 
applications. 

5.  Conclusions 

This  paper  investigates  on  the  thermoelectric  conversion  per¬ 
formance  of  SiC  ceramic  structures  on  supersonic-vehicle  nose  tip 
under  the  conditions  of  different  Mach  numbers.  The  numerical 
analysis  is  performed  with  and  without  respect  to  the  Thomson 
heat,  respectively. 

Results  show  that  there  is  better  thermoelectric  performance  at 
high  Mach  numbers.  With  the  increased  Mach  numbers,  both  the 
output  power  and  efficiency  become  higher  because  of  growing 
temperature  difference  and  heat  flux  attributed  to  the  dramatically 
increased  aero-heating.  The  maximum  output  power  reaches 
0.062  W  at  temperature  difference  of  670  K  with  M  4.5. 

When  the  Thomson  effect  is  considered,  the  thermoelectric 
performance  is  strongly  weakened.  There  is  over  50%  of  the  gener¬ 
ated  power  converting  to  Thomson  heat  due  to  the  large  tempera¬ 
ture  difference  at  each  Mach  number  case.  Notwithstanding,  the 


Current  [A] 

Fig.  13.  Output  power  and  energy  efficiency  versus  current  with  respect  to  Thomson 
effect  at  different  Mach  numbers. 


thermoelectric  conversion  performance  is  continuously  improved 
with  the  increase  of  Mach  number.  Moreover,  the  Thomson  heat  is 
possible  to  be  eliminated  by  the  optimization  of  Seebeck  coefficient, 
which  is  a  new  research  area  of  nanostructured  materials  and 
nanocomposites.  Therefore,  the  thermoelectric  technology  used  in 
aerodynamic  heat  recovery  is  very  promising. 

This  research  suggests  a  novel  application  of  thermoelectric 
uniform  structures  acceptable  for  harvesting  aerodynamic  heat 
energy  on  supersonic  vehicles,  and  also  supplies  fundamental  data 
for  limiting  the  aero-heating  by  thermoelectric  technology. 
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